Effective temperatures and activated dynamics for a two-dimensional air-driven 
granular system on two approaches to jamming 
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We present experiments on several distinct effective temperatures in a granular system at a se- 
quence of increasing packing densities and at a sequence of decreasing driving rates. This includes 
single-grain measurements based on the mechanical energies of both the grains and an embedded 
oscillator, as well as a collective measurement based on the Einstein relation between diffusivity 
and mobility, which all probe different time scales. Remarkably, all effective temperatures agree. 
Furthermore, mobility data along the two trajectories collapse when plotted vs effective temperature 
and exhibit an Arrhenius form with the same energy barrier as the microscopic relaxation time. 
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A critical challenge for the next decade is to under- 
stand nonequilibrium behavior such as commonalities in 
the jamming of glassy liquids, colloidal suspensions, and 
granular media [l|, [2| • The concept of effective temper- 
ature, defined through the relation between the fluctu- 
ations and the response of a nonequilibrium system to 
small perturbation Q, is an important unifying principle. 
Away from thermal equilibrium, behavior could depend 
on history and driving, so one would not expect a unique 
effective temperature [4]. Nevertheless, nearly ten dif- 
ferent definitions yield a common low-frequency value in 
simulations of steadily driven Lennard- Jones systems [H[ , 
foams [fij], and granular media [7]. Furthermore, acti- 
vated dynamics according to effective temperature have 
been incorporated into the theories of soft glassy rheol- 
ogy [g[ and shear-transformation zones Hi, a nd have also 
been demonstrated in simulation [l(J flJlll^ . 

Despite these advances, the utility of effective tem- 
peratures in real-world systems is not as clear. Viola- 
tions of the fluctuation-dissipation relation have been re- 
ported f or g lass-forming materials [H] and for granular 
systems [14| , however no more than one effective temper- 
ature has ever been measured and compared for a single 
sample. For aging colloidal glasses, there is even contro- 
versy about whether the Einstein relation between dif- 
fusivity and mobility gives an effective temperature any 
different from the bath temperature [151 ]. Furthermore, 
activated dynamics based on effective temperature have 
never been demonstrated in experiment. 

In this paper we report on a quasi-two dimensional 
granular system of macroscopic spheres, which roll with- 
out slipping on a horizontal plane. Uniform steady-state 
motion is excited by air, blown upward through a per- 
forated mesh on which the balls roll. Turbulent wakes 
are the source of both a rapidly-varying random driv- 
ing force and of a velocity-dependent drag force on each 
sphere [l6[; wake- wake interactions also cause a repul- 
sive force between neighboring spheres [17]. This system 
exhibits such hallmark features as a growing plateau in 
mean-squared displacement [18| and a growing dynami- 
cal correlation length [H, |2Q] on approach to jamming. 
Since the grains are macroscopic, and since the dynam- 
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FIG. 1: (Color online) (a) Distribution of kinetic energies 
K s and Kb, for small and big balls, respectively, and distri- 
bution of kinetic energy K t and potential energy U± for the 
weighted-ball thermometer, at air speed 670 cm/s and area 
density 0.43. Inset: scatter plot of instantaneous kinetic and 
potential energy for the thermometer, (b) Distribution of 
thermometer total mechanical energy. Inset: the thermome- 
ter, as seen from above, is a weighted hollow sphere whose 
position degrees of freedom are measured by the coordinates 
r of its geometric center and of a surface spot 5r marked ra- 
dially opposite the center of mass. The solid black curves in 
both (a,b) are the statistical mechanics predictions based on 
total average mechanical energy Tt = (Kt + Ut)/2. 



ics are slow, the ease of measurement and manipulation 
rivals that of computer simulations. Here we take advan- 
tage of this feature, both for comparing different effective 
temperatures as well as for evaluating their meaning in 
terms of activated dynamics. 

The experimental system consists of 50:50 bidisperse 
hollow polypropylene balls confined to a 12" diameter 
region; the diameters and masses are {2.54 cm, 2.2 g} 
and {2.86 cm, 3.0 g}. The balls are imaged from above 
at 30 Hz and 15 pixels/cm. Position and velocity data are 
extracted using custom Lab VIEW programs to respective 
accuracies of 0.001 cm and 0.03 cm/s. 
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The simplest effective temperature is based on kinetic 
energy, which has translational and rotational contribu- 
tions since the balls roll without slipping. Example re- 
sults for the kinetic energy distributions of both size balls 
are plotted in Fig. [Ha), for air speed 670 cm/s and ball 
area density 0.43. The two distributions are similar and 
roughly exponential, as in a two-dimensional thermal sys- 
tem. The average kinetic energy thus defines an effective 
"granular" temperature T g . At fixed airspeed, T g goes 
linearly to zero as the packing fraction increases toward 
0.78, similar to Fig. 9 of Ref. [18|. 

The second effective temperature is inspired by the- 
oretical consideration of a harmonically tethered test 
grain 0, [TO, Hl| . This has been realized experimentally 
by optical trapping of a probe particle in an aging col- 
loidal glass [22j |. Here a similar effect is achieved with a 
2.86 cm hollow plastic sphere that is partially filled with 
glue, so that the total mass is 3.8 g and the center-of-mass 
is 8.8 mm below the geometrical center. When tilted, 
the weighted-sphere "thermometer" stores gravitational 
potential energy and experiences a restoring force like a 
pendulum. The natural oscillation frequency is 2 Hz; the 
Q factor is about ten; the r.m.s. tilt is less than 20°, so 
oscillations are harmonic. A diagram of the thermome- 
ter is shown in Fig. [Jjb). The geometrical center of the 
thermometer and the location of the apex are tracked 
by video. From these quantities we compute the total 
kinetic energy K tl based on translational and rotational 
speeds plus mass and moment of inertia, as well as the 
total gravitational potential energy Ut, based on the rise 
in center-of-mass. 

Example results for the thermometer energies are plot- 
ted in Fig. [H The scatter plot of K t vs U t in the top 
inset reveals no correlations. The probability distribu- 
tions of K t and Ut in the top plot are nearly indistin- 
guishable, and appear to be exponential in reasonable 
agreement with the bath ball kinetic energy distribu- 
tions under identical conditions. This is characteristic 
of an object in thermal equilibrium with a heat bath, 
where the kinetic and potential degrees of freedom are 
independently populated and where the energies are ex- 
ponentially distributed according to a Boltzmann factor. 
This defines an effective temperature T t = (K t + U t )/2. 
The thermal nature of the thermometer motion is fur- 
ther illustrated in the bottom plot, Fig. []Jb), °f the 
distribution of total thermometer energy E = K t + U t . 
Since the system is two-dimensional, the density of states 
is proportional to E and the distribution should be 
P(E) = (E/T t 2 )exp(-E/T t ), which matches the data 
very well. The weighted sphere thus behaves like an ob- 
ject in thermal equilibrium, and T t is truly effective in its 
statistical mechanical meaning as a temperature. 

The third effective temperature is based on the Ein- 
stein relation that diffusivity D equals temperature times 
mobility \i. Since diffusion and flow involve grain- 
scale rearrangement of many neighboring beads, this 
probes long-time collective dynamics, as opposed to the 
short-time single-grain dynamics probed by T g and the 
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FIG. 2: (Color online) Mobility vs tilt frequency for 1" di- 
ameter test spheres at air speed 760 cm/s and area density 
0.57. Symbol types distinguish test spheres and tilt ampli- 
tudes; solid blue square: teflon, 18.45 g, 0.52°; hollow (solid) 
red diamond: acrylic, 10.12 g, 2.35° (0.52°); hollow (solid) 
green triangle: wood, 6.4 g, 0.52° (0.32°). Inset: instanta- 
neous velocity vs instantaneous force for the wood test sphere 
at 1.67 x 10" 3 Hz and tilt amplitude 0.52°. Mobility is the 
slope of the best fit to v oc F (black line). 



intermediate-time dynamics probed by T t . Here we mea- 
sure D from the linear growth of the mean-squared dis- 
placement over at least one decade in time. We measure 
\i for several 2.54 cm solid spheres with different masses, 
m, by tracking their motion while the entire bed is rocked 
sinusoidally in time. When the system is at angle away 
from horizontal, a test sphere experiences a force F down 
the plane given by mg sin 0, minus the buoyancy due to 
density difference with bath balls, plus a contribution set 
by the in-plane acceleration due to the location of rota- 
tion axis at the bottom of the windbox about four feet 
below the balls. Therefore, parallel to the tilt direction, a 
test sphere acquires an average drift speed of v = F\i that 
is superposed on its random thermal motion. To extract 
mobility, we therefore make a scatter plot of parallel in- 
stantaneous speed vs instantaneous force at all times and 
then we fit for the proportionality constant, as in the in- 
set of Fig. El Mobility values for several test spheres are 
plotted vs tilt frequency in the main plot. If the stochas- 
tic motion of the bath balls is unperturbed, the mobility 
results are independent of test mass, tilt frequency, and 
tilt amplitude; this demonstrates linearity of response. 
Thus we reliably measure both mobility and diffusivity, 
from which we compute an effective "Einstein" temper- 
ature T e = D/fi. 

Results for the three effective temperatures are com- 
pared in Fig. [3] as the system is brought closer to jamming 
along two different trajectories. For increasing area frac- 
tion at fixed airspeed 670 cm/s, and also for decreasing 
airspeed at a fixed area fraction 0.57, the various effec- 
tive temperatures are found to decrease by roughly one 
order of magnitude. However, the more crucial feature 
of Fig. [3] is that the values of all three types of effec- 
tive temperature are in agreement to within uncertain- 
ties estimated from the imaging resolution and fitting 
accuracies. Close to jamming the Einstein temperature 
systematically drops below the others, but this may be 
because the growing sub-diffusive plateau causes an un- 



3 



10 



a. U = 670cm/s - 



▼ T t 
■ T p 




0.46 0.52 0.58 0.64 0.70 

area fraction 



900 800 700 600 500 400 

airspeed [cm/s] 



FIG. 3: (Color online) Effective temperatures vs (a) area frac- 
tion at airspeed 670 cm/s, and (b) vs driving airspeed at area 
fraction 0.57. The granular temperature T g is based on the 
ball kinetic energies; the thermometer temperature Tt is based 
on the mechanical energy of weighted-sphere oscillators with 
two different masses, 3.83 g for upward triangles and 6.82 g 
for downward triangles; the Einstein temperature T e is the 
ratio of diffusivity to mobility. The solid green curve is the 
weighted average of all these measures. Measurement uncer- 
tainties for T g and Tt are comparable to symbol size, while 
the uncertainties in T e are indicated by error bars. 



derestimate of D. But overall, the good agreement allows 
us to average together {T gi T tl T e } with weights given by 
the uncertainties. The result is a single effective temper- 
ature, T e ff, displayed by solid green curves in Fig. [3j 

We now examine the mobility, not as a function of 
the actual experimental control parameters of area frac- 
tion and airspeed, but rather as a function of the ef- 
fective temperature. As usual for glass-forming liquids, 
we make a semi- logarithmic plot of l//i vs 1/T e ff in 
Fig. UK a). The values of T e ff/D, also included, differ 
slightly from because the Einstein temperature is 
not identical to T e ff . Here the temperature axis is scaled 
by ball weight times diameter, mgD^ the mobility axis 
is scaled by fi = r c /m, where r c is the cage-recognition 
time when grains cross from ballistic to subdiffusive mo- 
tion, as defined in Fig. 11 of Ref. [Hj]. Unexpectedly, 
we find that the data along the two different trajecto- 
ries collapse when plotted according to effective temper- 
ature. Scaled T g /D results for larger collections of steel 
balls [H, Hfll also collapse nicely, though with more scat- 
ter. The "thermodynamic" state of air-driven grains, and 
their distance from jamming, are entirely captured by the 
value of T e fi/ (mgDb) without regard to the actual control 
parameters of airspeed and density. Next, remarkably, we 
find that 1/ fi is exponential in E/T e R and hence is Ar- 
rhenius, just like in strong glass- forming liquids where re- 
laxation is activated with a single energy barrier E. Here 
the value of mgDb/E, indicated by the vertical gray line 
in Fig. [H is set by air-mediated repulsion of a grain by 
its surrounding cage of neighbors. Note \i « \i Q when 
T e ff w E, supporting the Arrhenius interpretation. 

While mobility can be thought of as a relaxation rate, 
it is instructive to compare with more direct microscopic 
time and length scales. In particular, the proximity to 
jamming may also be gauged from the mean-square dis- 
placement by the separation of the cage-recognition time 
r c and the rearrangement time r r , when grains cross 
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FIG. 4: (Color online) (a) Inverse mobility scaled by \i — 
Tc/m, (b) dimensionless extent of subdiffusive plateau in 
mean-squared displacement, and (c) average number of grains 
in kinetic heterogeneities, all vs inverse effective temperature 
scaled by ball weight times diameter. Open and closed sym- 
bols are for the polypropylene balls used in Figs. 1-3, along 
two trajectories as labeled. The x symbols are based on ear- 
lier data |l8[ for a bidisperse mixture of 0.635 and 0.873 cm 
diameter steel balls, and the + symbols are based on earlier 
data [2(| for a bidisperse mixture of 0.318 and 0.397 cm diam- 
eter steel balls; for both, jamming is approached by change 
in packing fraction at constant airspeed. In (a) the diamond, 
x, and + symbols represent T e ^/D. In (a-c) the dashed lines 
are Arrhenius fits, oc exp( J E/T e ff), all with the same energy 
barrier indicated by the vertical gray line. The solid curves 
are power-law fits fi /l^ oc r r /r c oc (mgDb/T e ^) 2 ' 0±0 ' 5 and 
n* oc (mgD b /T eE )°- 7±0 - 2 



from subdiffusive to diffusive motion. The rearrange- 
ment dynamics are spatially heterogeneous, punctuated 
by sudden motion of a string-like cluster of neighbor- 
ing grains [23[ whose average number n* may be de- 
duced from four-point correlation functions [20[. Previ- 
ously we measured both r r /r c [H, Hfll and n* [20j] vs 
packing fraction at fixed airspeed for larger collections 
of smaller steel balls. Now we compute these quanti- 
ties for all systems, and plot the results vs mgD^/T^ 
in Figs. H|(c,d). As m 9Db/T e & increases on approach to 
jamming, the extent r r /r c of sub-diffusive motion and 
the size n* of heterogeneities both collapse for the differ- 
ent grain sizes and also appear to grow without bound. 
Arrhenius fits with the same energy barrier as for mo- 
bility are satisfactory for r r /r c , but not for n*, both 
as expected. Power-law fits are also satisfactory, as 
shown. The latter give r r /r c ~ (mgDb/T^ ) 2 - 0±0 - 5 and 
n* ~ (m^^/Teff) ^ 0,2 ; these exponents are consistent 
with simulation of a Lennard- Jones liquid [23 |. Intrigu- 
ingly, the size exponent is also consistent with simulations 
of athermal systems vs packing density [25^ . 

In summary, we have presented the first experimental 
evidence that distinct effective temperatures can agree. 
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This holds across a range of conditions, and includes 
measures based on the kinetic energies of two sized bath 
spheres at short times, the kinetic and potential energies 
of two weighted-ball oscillators at intermediate times, 
and the ratio of diffusivity to mobility at long times. 
For a nonequilibrium system there is no a priori rea- 
son for these seven quantities to agree. Furthermore we 
have demonstrated that T e ff acts as the sole state vari- 
able for air-driven grains and that relaxation times are 
Arrhenius in T e fj. Even more than the measures of struc- 
ture and dynamics in [l8[ , and the measures of spatially- 
heterogeneous dynamics in pjl [20[ , this establishes air- 
fluidized beads as a faithful model for the glass transi- 



tion. Besides extending the universality of the jamming 
and effective temperature concepts, this is also important 
because key microscopic quantities such as r r /r c and n* 
are not as fully accessible in other experimental systems. 
One interesting line of future research would be to im- 
pose shear, or an unlikely initial configuration, so that 
the system is away from "equilibrium" , and to study dif- 
ferences in the various effective temperatures. An even 
broader line would be to explore when and why driven 
systems exhibit equilibrium-like behavior [26]. 
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